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ABSTRACT  This  paper  analyzes  the  capability  of  low  orbit 
satellite  borne  interception  devices  to  carry  out  orbital 
maneuvers  to  intercept  high  orbit  satellites.  In  interception 
system  analysis,  threat  areas,  defense*  areas,  and  intercept  areas 
are  three  basic  concepts.  This  article  gives  methods  for  their 
modeling  as  well  as  solution.  In  conjunction  with  this,  use  is 
made  of  digital  simulation,  giving  standard  threat  areas,  defense 
areas,  and  intercept  areas  corresponding  to  different  intercept 
time  periods. 


KEY  TERMS  Satellite  orbit.  Space  flight  interception 
device.  Intercept  trajectory 
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I .  THE  CONCEPTS  OF  THREAT  AREAS ,  DEFENSE  AREAS ,  AND 
INTERCEPT  AREAS 
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Interception  of  high  orbit  satellites  by  low  orbit 
satellite  borne  interception  devices  is  a  "space  based-space 
based"  intercept  system.  Intercept  systems  are  composed  of 
satellite  platforms  moving  along  respective  orbits  as  well 
as  interception  devices  I  capable  of  orbital  maneuver  and 
carried  on  the  platforms.  Intercept  targets  are  high  orbit 
satellites  SS.  The  intercept  debris  path  is  as  shown  in 
Fig.l.  At  the  instant  when  t*=0,  SS  moves  along  the  orbit 
to  point  S.  At  the  instant  t  (illegible)  arrives  at  point 
S.  At  instant  T(illegible)  the  satellite  platform  moves  to 
point  B.  At  this  time,  the  satellite  platform  releases  I. 

In  accordance  with  guidance  methods,  I  carries  out  orbital 
maneuver  along  a  certain  intercept  orbit  and  also  arrives  at 
point  S(illegible)  at  just  the  instant  T(illegible^  The 
point  S*  is  called  the  predicted  hit  point.  At  this  point, 

I  collides  with  SS,  destroying  the  target. 

important  concepts  in  intercept  system  analysis  are 
threat  areas,  defense  areas,  as  well  as  intercept  areas. 
Definitions  are  as  follows* 


If  one  assmieB  that  the  instant  when  intercept  begins 
is  r(illegible)-0,  at  this  time,  SS  is  located  at  point  S, 
and  the  satellite  platform  is  located  at  point  B.  Due  to 
the  fact  that  I  orbital  intersection  maneuver  capabilities 
are  limited  by  the  amount  of  fuel  carried  as  well  as  other 
restrictions,  speaking  in  terms  of  point  S,  therefore,  there 
exists  a  zone  where,  if  only  the  satellite  platform  is 


*  Numbers  in  margins  indicate  foreign  pagination. 
Commas  in  n\imbers  indicate  decimals. 
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positioned  within  this  zone  at  instant  t (illegible )»0,  and, 
in  conjunction  with  this,  I  is  released  at  this  instant, 
going  into  orbital  change  maneuvers,  it  flies  along  a 
collision  course.  Only  then  is  it  possible  to  make  I 
collide  with  SS  at  the  instant  r*  .  This  zone  is  called  the 
T*  threat  zone  of  SS.  Obviously,  satellite  platforms 
outside  T (illegible)  threat  zones  do  not  possess  the 
capability  of  killing  or  damaging  SS  located  at  position  S 
at  instant  t*=0.  The  greater  the  number  of  satellite 
platforms  is  within  the  t*  threat  area,  the  greater,  then, 
is  also  the  number  of  I  possessing  the  capability  of  killing 
or  damaging  SS  at  time  r*.  As  a  result,  the  greater  is  also 
the  threat  posed  to  SS.  Because  of  this,  speaking  in  terms 
of  satellite  platforms,  threat  areas  determine  whether  or 
not  the  I  they  carry  possess  the  capability  to  kill  or 
damage  SS.  Speaking  in  terms  of  SS,  threat  areas  determine 
the  magnitude  of  the  threats  they  pose  as  well  as  the 
locations  on  their  own  orbits  giving  rise  to  collisions 
(determined  by  t (illegible) ) . 

In  the  same  way,  with  regard  to  satellite  platforms  at 
instant  t*=0  located  on  point  B,  there  also  exists  a  zone. 

Only  when  SS  enters  this  zone  at  instant  r*  is  it  then 
possible  to  be  killed  or  damaged  by  platforms.  This  zone  is 
called  the  platform  r*  defense  area,  t*  defense  areas 
determine  the  size  of  platform  operational  ranges,  that  is, 
satellite  platforms — after  given  intercept  time  periods 
T-*— are  only  capable  of  killing  or  damaging  SS  that  enter 
the  areas  in  question.  However,  intercept  areas,  by 
contrast,  are  defined  as  the  collection  of  all  possible 
intercept  collision  points  (that  is,  predicted  hit  points) 
for  a  certain  platform  and  intercept  time  period  t*.  From  a 
different  angle,  they  reflect  the  size  of  platform  operation 
ranges . 
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T*  threat  areas ,  r*  defense  areas,  and  t*  intercept  areas — 
besides  being  related  to  intercept  time  periods  t* — are  also 
related  to  the  locations  of  platforms  and  SS  at  instant  r*,  fuel 
by  interception  devices,  as  well  as  guidance  methods*  Xn 
this  article,  option  is  made  for  the  use  of  speed  gain  guidance 
to  act  as  the  guidance  method  for  I. 

This  article  will  give  precise  methods  for  the  modeling  and 
solution  of  these  three  types  of  zones.  In  conjunction  with 
this — going  through  digital  simulation — it  will  present  classical 
T*  threat  areas,  t*  defense  areas,  and  t*  intercept  areas. 


II.  RESEARCH  ASSUMPTIONS,  COORDINATE  SYSTEMS,  AND  SATELLITE 
PLATFORM-I-SS  KINEMATICS  MODELS 

( I )  As  sumptions 

The  entire  interception  process  associated  with  low  orbit 
satellite  borne  interception  devices  intercepting  high  orbit 
satellites  occurs  outside  dense  atmospheric  layers.  In  research, 
the  assumptions  made  are  as  follows: 

1.  Spacecraft  fly  outside  the  atmosphere.  Satellite 
platforms  and  SS  only  undergo  the  effects  of  the  earth's  gravity. 
During  flight  processes,  orbital  maneuvers  are  not  carried  out. 
Besides  undergoing  the  effects  of  the  earth's  gravity,  I  also 
undergo  the  effects  of  thrust.  Thrust  uses  the  form  of  impulses 
at  points  of  orbital  change  to  influence  I,  making  I 
instantaneously  achieve  required  velocity  gains. 

2.  The  globe  is  a  sphere  with  radius  R( illegible)  — 
.6371110m.  The  gravitational  field  of  the  earth  is  an  inverse 

square  of  the  distance  force  field.  The  gravitational  constant 
of  the  earth  is  ^=3.986005xl014m3/s2. 
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3.  Disregard  the  effects  of  global  rotation  on  spacecraft 
movements . 

4.  Satellite  platforms  are  circular  polar  orbits  with 
altitudes  of  H=640km.  High  orbit  satellites  are  stationary 
equatorial  satellites  with  orbital  altitudes  of 
H(illegible)=35787km.  Maximum  velocity  impulses  which 
interception  devices  are  capable  of  supplying  themselves  are 
■Vmax=6421m/s . 


(II)  Coordinate  Systems  and  Coordinate  Transforms 


1.  In  the  geocentric  inertial  coordinate  system  OeXYZ,  Oe 
is  congruent  with  the  center  of  the  earth.  The  X  axis  in  the 
equatorial  plane  points  toward  the  vernal  equinox.  The  Z  axis 
and  the  direction  of  autorotation  of  the  earth  are  the  same.  The 
Y  axis  is  determined  by  a  right  hand  rule,  as  shown  in  Fig. 2. 


2.  As  far  as  satellite  platform  orbit  coordinate  system  Oe- 
XpYpZp  is  concerned,  Xp  is  the  direction  of  continuation  of  Oe 
and  the  instant  t*=0.  Yp  lies  in  the  satellite  platform  orbital 
plane,  points  in  the  direction  of  flight,  and  is  positive.  Zp  is 
determined  by  a  right  hand  rule  as  shown  in  Fig. 2. 

3.  With  regard  to  interception  orbit  coordinate  system 
OeXIYIZI,  XI  and  Xp  are  congruent.  YI  lies  in  the  plane  of 
interception  orbit,  points  toward  the  predicted  target  hit  point 
S*  ,  and  is  positive.  ZI  is  determined  by  a  right  hand  rule  as 
shown  in  Fig. 2.  In  order  to  facilitate  the  researching  of 
problems  one  may  as  well  assume  that  the  position  under  the  stars 
of  the  vernal  equinox  is  just  congruent  with  S*  .  Then,  at  the 
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instant  t*=0,  SS  should  be  positioned  at  the  point  S  on  the 
eastern  equatorial  radius  -nasT*.  In  this. 


is  the  orbital  angular  velocity  of  SS. 


4.  Relationships  Between  Oe-XYZ  and  Oe-XpYpZp. 


Taking  Oe-XYZ  and  first  turning  it  around  the  Z  axis  an 
angle  £J,  then  turning  it  around  the  new  Y  axis  -uo,  and  finally 
turning  90*  around  the  new  X  axis,  it  is  then  possible  to  obtain 
Oe-XpYpZp.  Coordinate  transform  relationships  then  are: 


(1) 


cosQcosz/^ 

—  cosfisinMjj 

sinQ  1 

sinQcosu,! 

—  sinOsint/,, 

—  cosO 

sinw„ 

.  cost/ II 

0 

In  this,  n  and  uo  are  the  longitude  and  latitude  associated  with 
I  at  instant  t*=0.  /41 

5.  Relationships  Between  Oe-XIYIZI  and  Oe-XpYpZp. 

The  included  angle  between  intercept  orbit  planes  and 
satellite  platform  orbit  planes  are  a.  Turning  XpYpZp  around  Xp 
by  a,  one  obtains  Oe-XIYIZI.  The  coordinate  transform 
relationships  then  are: 
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In  thlS/  oc  is  cnpabie  of  being  solved  for  front 
the  relationships  set  out  below: 

-90  °-<u„< 0  °  ,  tg£2  =  tg3tsin(-M„);  / 

270  °  <w<,<“90  °  ,  ctga  =  ctgficos(-M„-90  ^ );  ■ 


For 

For 


(III)  Positions  and  Speeds  of  Satellite  Platforms  at 
Instant  r* 


In  the  system  Oe-XpYpZp,  one  has: 


(4) 


v^=io,4JJjrTH),of 


(IV)  I  Kinetics  Models 


In  the  Oe-XIYIZI  system,  I  movement  differential  equations 

are: 
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(5) 


[  ^ I  ^  ^Tx 

\Y,~-„y,/r'+a„ 

lz,=  -pr,/i-’+<-„ 


In  this: 
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r  =  ^x]  +  Y]+z] 


aTX,  aTY,  and  aTZ  act  as  thrust  accelerations  used  on  unit 
masses.  f 

When  thrust  accelerations  are  0,  I  is  in  free  flight. 

At  the  instant  t*=0/  satellite  platforms  release  I.  Due  to 
the  fact  that  the  XI  and  Xp  axes  are  congruent,  as  a  result,  the 
initial  I  location  vector  and  the  initial  satellite 

platform  vector  are  the  same.  F,  =  =  (/?^  +  y/,0.0) 

If  one  adopts  thrust  impulse  assumptions  that  I,  at  the  instant 
T*=0,  achieves  velocity  increase  amounts  AK  and  record  is 
made  of  AF  in  ttie  system  Oe-XI  YI  ZI  expressed  as  AF 

=  ^  >  then,  the  initial  velocity  in 

system  Oe  -  XI  YI  ZI  can  be  expressed  as: 
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(VI)  SS  Kinematic  Models 


Assuming  that,  at  the  instant  t*=0,  SS  is  located  at 
position  S  and  flies  at  the  instant  t*  just  to  position 
S(illegible) ,  then,  the  geocentric  angle  spread  between  S  and 
S( illegible)  is  (illegible).  S( illegible)  is  the  predicted  hit 
point.  At  the  instant  r (illegible) ,  the  SS  position  vector  in  Oe 
-XYZ  is  (  A’.t//..().())^  ,  then,  in  the  system  Oe- 

XIY( illegible )ZI,  expressed  as: 


(7)  '(A,  -r 

III.  IMPULSE  VELOCITY  GAIN  GUIDANCE 

In  the  intercept  orbit  coordinate  system  Oe  -  XI  YI  ZI,  the 
I  position  vector  is  The  SS  position  vector  is  /T. 
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Then,  delta  f  is  already  known,  as  shown  in  Fig.l.  At  instant 
T( illegible)  =0,  satellite  platforms  release  I.  Then,  I  possess 
the  platform  velocity  i’„  .  In  the  system  Oe  -XI  YI 

Z (illegible) ,  it  is  expressed  as  .  After  that,  by 

thrust  impulses  supplied  instantaneously  by  propulsion  systems  it 
carries  itself,  I  is  made  to  achieve  a  velocity  gain  of  AT 
Then,  I  changes  instantaneously  from  velocity 

'P‘,  .  ,  F  is  designated  as  the  necessary  velocity, 

that  is,  with  I  flying  along  the  intercept  path  and — at  the 
instant  t* — hitting  the  predetermined  hit  point  £*,  the  required 
velocity  at  instant  t*  =0.  After  solving  for  ,  one 

obtains  . 

-^/o)  )-  described  above). 

Due  to  F,  r  it  is  possible  to  surpass  the  second  cosmic 
velocity.  As  a  result,  intercept  orbits  can  be  a  certain  type  of 
circular  conic  curve — ellipse,  parabola,  or  h3^erbola. 

Below,  use  is  made  of  Herrick  methods  in  order  to  precisely 
determine  intercept  orbits  and  required  velocities  Vr  .  Herrick 
methods  first  surmise  that  intercept  orbit  semiparameters  Po 
associated  with  setting  out  from  point  B  at  instant  t*  and 
passing  through  point  S*  are — due  to  ^  ,  r,  ]  already 

known  quantities.  Therefore, 


(8) 


cosA/  =  /• 


//•/, 


/43 
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(9) 


r 


=  /V(l+ecos/J  /j=l,2 


One  then  has 


ecos4  =/>/r^ -1  n=l2 
es'mf ,  =  (ecos/,  cosA/-  ecos/,)sinA/ 
esin/,  =  (ecos/  |  -ecos/,cosA/)sinA/ 


Therefore,  it  is  possible  to  solve  for 


J\^  =  arctg(esin4  /  ecos/’^ )  «  =  1,2 


Due  to 
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n  =  1.2 


Therefore 


Make 


(10) 


sin£ .  — sin/,, 


•  sl—e'  /  P 


>cos£.,  =/•„(<' 

sin/;//  Mil/  ,  \  /  P 

n=\.2. 

[  cos/;  II  '■  ,  i ^ni>/  „ )  •'  P 


\  I 

*;  ellipse ! 

••i 

I 


I 

t'  >  ij  hyperbole ; 


one  has 


f  £  =  arciciMn/: ,  /  cos£ ,  )  /;=1.2  t- <  ^  ellipse 

•j  // ,  =  arcig/;(sin/;//„  /  cosliH^  )  ;i  =  1.2  e>  1  hyperbole 
if  /2)  h  =  1.2  t'=l  parabola 

II  •^^11 

a  =  P/{l-e') 


A/„=£  -esin£ 
=esin////^-//^_ 

=(^„+/=';/3)/2 


/I  =  1,2 

;;=  1.2 

«=  1.2 


e<t  ellipse 

^  ^  ^  h^erbole 
e =  1  parabola 
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Then 


(11) 


T* 

H  t’  =V 

T*  =TPV^{M  M 


g  <  i',  ellipse 
e  >  hyperbola 
e  =  parabola 


If  the  T*  solved  for  and  given  values  are  not  equal,  then, 
corrections  are  made  to  guessed  P  values,  and  the  operations 
described  above  are  repeated  until  the  difference  between  r* 
solved  for  and  given  values  satisfies  accuracy  requirements  and 

they  are  stopped. 

In  a  case  where  intercept  orbit  semiparameter  P  is  already 
known,  it  is  then  possible  to  solve  for  Vr  components  in  the 
system  Oe  ~  XI  YI  ZI  as: 


-X,  t' 

(f 

£  -±- 

-  l)cosA/-(p  - 1) 

V  = 

r 

*'• 

P 

sinA/ 

V/xP  /  f  j 

_ 

0 

(12) 
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Thus,  the  size  of  required  velocity  gains  can  be  precisely 
determined  from  the  equations  below: 


(13) 


AV=V  -B''v 

r  0 

AK=  |a7| 


Due  to  limitations  on  the  fuel  I  carries,  if  the  maximum 
velocity  impulse  I  propulsion  systems  are  capable  of  supplying  is 
delta  Vmax  ,  then  it  is  required  that 


(14) 


A  I'  <  A  F 

mii\ 


Besides  this,  attention  should  also  be  paid  to  the  fact  that 
I  is  only  suitable  as  an  interception  device  for  flight  outside 
the  atmosphere.  If  the  interception  orbit  during  the  process  of 
flying  from  B  to  S  has  a  minimum  altitude  from  the  surface  of  the 
earth  of  Hm  ,  then  it  is  required  that 


(15) 


>  lOOkm 
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IV.  METHODS  TO  PRECISELY  DETERMINE  T*  THREAT  AREAS,  T*  DEFENSE 
AREAS,  AND  T*  INTERCEPT  AREAS 


(I)  Precise  Determination  of  t*  Threat  Areas 
During  the  process  of  SS  intercept,  I  undergoes  restrictions 
associated  with  the  two  equations  (14)  and  (15).  As  a  result, 
there  are  only  I  carried  on  satellite  platforms  within  a  fixed 
range  which  then  possess  the  capability  to  kill  or  damage  SS  at 
instant  t*  .  This  is  the  threat  area  with  regard  to  SS.  Due  to 
satellite  platforms  making  circular  orbital  movements,  threat 
areas  are,  therefore,  a  section  of  spherical  surface  using  (Re+H) 
as  its  radius. 

At  instant  r*  =0,  the  angular  location  of  satellite ^platform 
B  relative  to  S  can  be  described  using  the  two  angles 
See  Fig. 2.  Using  the  center  of  the  earth  as  the  center  of  the 
sphere  and  (Re+H)  to  be  the  radius,  a  sphere  is  made.  The 
projection  on  it  of  the  north  pole  of  the  earth  acts  as  the  north 
pole  of  this  sphere.  SS  on  the  projection  is  S.  The  projection 
of  the  predicted  hit  point  is  S*.  A  meridian  line  is  made 
through  S.  Again,  through  B,  make  a  great  circle  perpendicular 
to  the  meridian  line  at  B’ .  Let  the  spherical  surface  angles  2 

BB'  =f5i/,and  SB'  =.Sq},  be  used  to  describe  the  location 

of  satellite  platform  B  relative  to  S.  Based  on  a  knowledge  of 
spherical  triangles,  it  is  not  difficult  to  obtain  the 
relationships  between  Su,  Sw  and  Q,  uo  ,  and  r*  . 

When  -90'<«„<0'’  ' 
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fi = arccos[cos«„cos(Q+/i.„T ' )] 
/BSB'  =90  ° +arcsin(sin«„/sin^) 


^4 
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^«,  =  arcsin(sinZBSB '  sin/J)  •  sing(Q) 
dll = — arccos(cos^  /  cosdoi) 


When  -270*  <m„<-90  * 


3(0  =  — arcsin[sin(£5+/j„T  *  )sin(— 90  °  )]  i 
5m =-90  *  -arccos[cos(— Mo-90  *  )/cos5co)] 


Using  5m.  to  be  the  transverse  coordinate  and  5a)  to  be 
the  longitudinal  coordinate  while  stipulating  north  of  the 

equator  as  positive  and  5a)  east  of  the  meridian  plane  to  be 
positive,  it  is  possible  to  solve  for  the  r*  threat  area 
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associated  with  platforms  in  direct  motion  is  as  shown  in 
Fig. 3 (a).  The  t*  threat  area  associated  with  platforms  in 
retrograde  motion  is  as  shown  in  Fig. 3(b). 


(II)  Precise  Determination  of  t*  Defense  Areas 


After  the  given  intercept  period  t*,  for  each  given  uo,  it 
is  possible  to  obtain  angular  distances  S  relative  to  satellite 
platform  orbit  surfaces  as  (fi+nssr*).  Because  £2  exists  within  a 
certain  range,  (Q+nssr*),  therefore,  also  has  an  interval.  It  is 
displayed  as  a  section  of  circular  arc  on  the  equator  (because  SS 
is  only  capable  of  moving  on  the  equator).  This  is  the  t* 
defense  area.  Defining  uo  north  of  the  equator  as  positive,  S 
east  of  the  platform  orbit  surface  as  positive,  and  taking 
(D+nssT*)  as  the  longitudinal  coordinate  and  Uo  as  the  transverse 
coordinate,  it  is  then  possible  to  obtain  r*  defense  areas 
associated  with  a  series  of  uo  as  shown  in  Fig. 3(c)  (direct 
motion  platform)  and  Fig. 3(d)  (retrograde  platform). 


(Ill)  Precise  Determination  of  r*  Intercept  Areas 


Intercept  areas  are  collections  of  all  possible  predicted 
hit  points.  On  the  equator  they  are  also  displayed  as  a  circular 
arc  section.  Positions  of  predicted  hit  points  relative  to  I  are 
described  using  C  and  uo.  Each  given  individual  uo  corresponds 
to  a  range  associated  with  a  certain  fi.  Taking  R  to  be  the 
longitudinal  coordinate  and  uo  to  be  the  transverse  coordinate. 
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it  is  then  possible  to  obtain  t*  intercept  areas  associated  with 
a  series  of  uo  as  shown  in  Fig. 3(e)  (direct  motion  platforms)  and 
Fig. 3(f)  (retrograde  motion  platforms). 


CONCLUSIONS 


1 .  As  far  as  interception  devices  which  are  carried  by  low 
altitude  satellite  platforms  with  orbital  altitudes  of  640km  and 
possessing  velocity  gains  of  roughly  6km/s  are  concerned,  they 
are  capable  of  intercepting  equatorial  stationary  satellites.  In 
accordance  with  current  technology  levels,  this  type  of  small 
interception  device  is  capable  of  realization.  As  a  result, 
equatorial  stationary  satellites  face  a  real  threat. 

2 .  There  is  a  relationship  between  the  size  of  threat 
areas,  defense  areas,  and  intercept  areas  and  intercept  time 
periods.  The  smaller  t*  is,  the  smaller  the  zones  then  are. 
Conclusions  clearly  show  that,  when  intercept  time  periods  are 
smaller  than  4x103  s,  interception  capability  still  exists. 
However,  within  this  time  period— with  the  help  of  advanced 
warning,  reconnaissance,  and  other  similar  means  SS  still  has 
adequate  time  to  adopt  a  number  of  effective  counter  interception 

measures . 


3.  The  methods  and  conclusions  supplied  by  this  article  can 
be  useful  in  the  analysis  of  space  warfare  systems. 


A 
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